Quantum synchronization and transresistance quantization in superconducting devices 
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We show theoretically the possibility of quantum synchronization of Josephson and Bloch oscilla- 
tions in a superconducting device. One needs an LC oscillator to achieve exponentially small rate 
of synchronization errors. The synchronization leads to quantization of transresistance similar to 
that in (Fractional) Quantum Hall Effect. 



One of the most interesting discoveries of XX cen- 
tury was the perfect (fractional) quantization of Hall 
transresistance in rather imperfect 2DEG semiconduct- 
ing samples [T]. The resistance as a function of electron 
density and magnetic field tends to be close to plateaus 
with values 
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n, to being integer numbers. The accuracy is so good as 
to enable numerous metrological applications [2, 3]. The 
physical explanation of the effect is the commensurabil- 
ity of electron density and density of the magnetic flux 
penetrating the sample, this taking place any time the 
ratio of numbers of elementary charges and flux quanta 
in the structure is a rational fraction 77/777. 

Quantum Hall samples are macroscopic involving in- 
finitely many degrees of freedom. Shortly after the dis- 
covery, Likharev and Zorin [4j hypothesized that sim- 
ilar resistance quantization may occur in a Josephson- 
junction superconducting device encompassing only few 
quantum degrees of freedom. They foresaw it as a result 
of quantum synchronization of Bloch [5] and Josephson 
[B] oscillations in two junctions. The Josephson frequency 
ljj = 2eVo /h is proportional to the average voltage drop- 
ping at one of the junctions while the Bloch frequency 
ljb = t^Iq/c is proportional to the average current in an- 
other junction. A synchronization condition of the two 
oscillations, nujj = mujB results in 

R= Y° = ^. (2) 
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The resistance quantum is modified in comparison with 
Eq. [T] manifesting the double charge 2e of Cooper pairs 
in superconductors. Unfortunately, the original device 
suggestion 11] does not work. The reason of the failure 
seems fundamental. The quantities to be synchronized, 
the charge and flux in the device are canonically con- 
jugated variables. Quantum mechanics forbids them to 
be simultaneously certain, and the synchronization is ex- 
pected to be destroyed by quantum fluctuations. 

A recent outburst of theoretical and experimental ac- 
tivities concerns quantum-coherent phase slips in thin 
nanowires [7j. On theoretical side, a concept of phase-slip 
(PS) junction has emerged |9]. Such junction is ex- 
actly dual to a common Josephson junction with respect 



to charge-flux conjugation. This inspired the proposals 
of novel superconducting devices [T0HT2]. Very recently, 
a PS qubit on InO nanowires has been realized [32] . Rel- 
evant experimental developments include observation of 
the predicted phenomena: phase-slips in Josephson junc- 
tion chains [HI HI], Bloch oscillations [16], and charge 
sensitivity [17] . 

In this Letter, we demonstrate that combining PS and 
Josephson junctions in a single device solves the prob- 
lem of quantum synchronization. A necessary element 
of the device appears to be an LC oscillator with high 
quality factor Q. With this, one can make the rate T of 
synchronization errors exponentially small, — lnT ~ Q 
thereby achieving exponential accuracy of the resistance 
quantization. Importantly, the device suggested can be 
also used as both voltage and current standard, thereby 
closing the metrological triangle [18 . 




FIG. 1. A PS (marked with diamond) and Josephson (cross) 
junction embedded into a general linear circuit. The circuit 
parts in dashed boxes generate Bloch and Josephson oscil- 
lations while the (frequency-dependent) resistors Z\ and Z2 
provide the coupling between the parts. The circuit is con- 
trolled with voltage and current sources Vb, h- The d. c. out- 
put voltage and current Vo,Io manifest the quantized tran- 
sresistance R—Vo /Io ■ 

To appreciate the difficulty of quantum synchroniza- 
tion, we consider first a PS and a Josephson junction 
embedded in a general linear circuit that can be rep- 
resented with four (frequency-dependent) resistors (Fig. 
[lj. The circuit parts in the dashed boxes represent the 
Bloch and Josephson oscillators. Let us first consider 
them separately by setting two coupling resistors Z\ t 2 to 
Z\ = 00, Z 2 = 0. The Josephson part is then a com- 
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mon (Til) Josephson current-biased junction shunted by 
the conductor Gj. If the bias current exceeds the criti- 
cal one, lb > Ic = 2eEj/H, the circuit produces voltage 
oscillations with frequency loj — 2e \° = ^f- y/I 2 — 1'^,, 
Vo being the time-averaged voltage across the junction. 
The energy accumulated in the oscillation is of the order 
of Josephson energy Ej. To have a well-defined semi- 
classical oscillation, we shall require that the energy ac- 
cumulated by far exceeds the quantum frequency scale 
fiu j. The latter can be regarded as an effective noise 
temperature TJ characterizing the quantum fluctuations 
in the circuit (we neglect the thermal fluctuations as- 
suming sufficiently small temperature). The condition 
Ej >• Tj amounts to Gj >• e 2 /H, the conductance must 
be high at quantum scale. 

The Bloch oscillator is understood with using the du- 
ality transformation between the phase and charge [S]. 
Upon such a transformation, the Josephson junction is 
replaced by a PS junction, the current bias by the voltage 
bias, and the parallel conductor becomes a series resistor 
Rs- Bloch oscillations occur provided the bias voltage 
exceeds the critical voltage of the junction, Vb > Vc = 
irEs/e. Their frequency uo B = ^ = ^yJV b 2 - V 2 , is 
related to Io, the average current in the junction. To 
have a well-defined semiclassical oscillation, we shall re- 
quire that the energy accumulated ~ Eg by far exceeds 
the effective noise temperature T B ~ Jtlob- This gives 
Rs 3> fi/e 2 : for a PS junction, it is the resistance that 
must be high at quantum scale. 

Let us now couple the circuits. The main effect of 
the coupling is the transfer of oscillating voltage/current 
from Josephson/Bloch to Bloch/ Josephson part, whereby 
the voltage/current is multiplied with the amplification 
coefficient K(lu) = Z 2 j[Z\ + Z 2 ). Besides, the effective 
resistance/conductance of Bloch/Josephson part is mod- 
ified, SR S = Z 2 Z 1 /(Z 2 + Zi), 5G.j = 1/(Z 2 + Z x ). In 
order to preserve well-defined oscillations, we require this 
modification to be small, SRs <C Rs, 5Gj <C Gj. 

We estimate the energy scale E sp associated with the 
coupling and synchronization of the oscillations as a 
product of oscillating voltage and current in each de- 
vice times oscillation period, assuming lob — ujj — uj, 
E sp ~ IqK{lo)Vq /uj. It is important to recall that 
the oscillating quantities are fundamentally related to 
frequency, Iq — eu>, Vo — homega/e. With this, 
E cp ~ KHuj. A generic estimation for K is K < 1. 
Indeed, for real impedances Z\ >2 K < 1. In this case 
Ecp ^5 

T B j and the envisaged synchronization in a gen- 
eral circuit is destroyed by quantum fluctuations. 

To overcome this, we need large K. An active amplify- 
ing circuit could provide this but brings extra noise that 
increases the fluctuations. The main idea of this Letter 
is to use a passive amplifying circuit, an LC-oscillator, 
replacing Z\ with a capacitor C and Z 2 with an inductor 
L (Fig. 2). With this, K(uS) ^> 1 near the resonant fre- 
quency £1 = (LC) 1 / 2 . Assuming that a small real part of 



Z 2 gives rise to a finite quality factor Q of the oscillator, 
K = (2(co/fl — l)+iQ)~ 1 at ui « f2. The maximum value 
of K is thus limited by Q, leading to E cp ~ Qhuj T} B . 
We expect the synchronization errors to be related to ac- 
tivation over this energy barrier and thus to occur at ex- 
ponentially small rate ~ exp(— E cp /T*) ~ exp(— aQ), a 
being a coefficient of the order of 1. We stress and prove 
further that the synchronization takes place in a rather 
broad interval of frequencies near f2: the Josephson and 
Bloch oscillations are thus synchronized with each other 
rather than with the LC oscillations. 
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FIG. 2. Coupling the Josephson and Bloch parts with an LC 
oscillator results in K ^> 1 close to the resonant frequency Q 
and thus enables quantum synchronization. 



The effective quality factor in our circuit is in fact lim- 
ited by dissipation in Rs,Gj. The conditions of non- 
obtrusive coupling SGj -C Gj,5Rg -C Rs imply that 
Q -C mm(GjZo, Rs/zq). In fact, the corresponding 
equality estimates the maximum effective quality factor 
Q~ = 1/GjZq + zq/Rs- A simple optimization of Q m is 
to chose the oscillator impedance zq = yj Rs/Gj, so that 
Qm = V RsGj/2. 

Synchronization persists in a finite interval of frequen- 
cies wb(H), ujj(Ib) near the line where those satisfy 
a given fractional ratio lub/ujj = n/m. To estimate 
the width of the interval, we compare E cp with an en- 
ergy scale characterizing the frequency deviation, which 
is either (Aujb/^b)Es or (Aujj/ujj)Es, this leads to 
(Auj b /ub) * K/(R s e 2 /h), (Auj/loj) ~ K/(Gjh/e 2 ) 
We see that for the limiting Q and at frequencies close to 
n the width of these intervals may become comparable 
with the frequency itself, Auij/ujj, Alu b /ub — L 

In the remainder of the Letter, we support these qual- 
itative estimations with quantitative illustrations. 

The adequate quantum description of the circuit in- 
volves two variables: superconducting phase drop at the 
Josephson junction 4> and dimensionless charge q = - Q 
flown in the PS junction. It is obtained in the framework 
of Keldysh action formalism 20 where variables are dou- 
bled <j) — > (t) , q —> q (t) corresponding to two parts of 
the Keldysh contour. It is convenient to use "classical" 
and "quantum" variables defined as 20, <f>d — ((f> + ±0"), 
2?,9d = {q + ± q~)- The total Keldysh action 



S — S B + Sj + S c 
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is contributed by the Bloch and Josephson parts 
Sj = f dt \ 2Ejsin<t>sm-^ - ^<^ + #d^ 



Sjj = I ill j 2E S sing sin y - ~^~ qd + qqd ~ 



the coupling part 

S C v — I ~ — 
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and the noise part that is quadratic in q d , 4> d and satisfies 
fluctuation-dissipation theorem (see [21] for concrete ex- 
pressions). The resulting action is non-local in time and 
therefore cannot be treated exactly. 

To start with, we study the resulting saddle-point clas- 
sical equations [21] neglecting the noise. This approxima- 
tion gives a good estimation of the positions and widths 
of the synchronization domains while disregarding round- 
ing of large and vanishing of small domains. Typical re- 
sults are presented in Fig. 3. For this plot, we made 
(mostly for esthetic reasons) a symmetric choice of pa- 
rameters Es = Ej, Gjh-K/Ae 2 = e 2 Rs/irfi, so that out- 
put current and voltage, and correspondingly the oscil- 
lation frequencies are symmetric in the plane of Vf, and 
lb, uj B {h/Ic,V b /Vc) = uj(Vb/Vc,h/Ic)- In average, 
these frequencies are close to those of uncoupled oscil- 
lators, u)B(Vb),u>j(Ib), the deviations are mostly due to 
synchronization. We observe the domains corresponding 
to the fractions n/m. They are centred at the curves 
where mws(H) = nu>j(h)- The widest domain is that 
with n — 1, m = 1 and is centered at the diagonal. The 
domains with higher n, m are increasingly more narrow, 
as it is also the case in QHE. The parameters are cho- 
sen such that the resonant frequency £1 is achieved at 
h/Ic — Vb/Vc — v2, were the domains are widest. 
Rs = I0irh/e 2 and the oscillator impedance is optimized, 
z o = y/Rs/Gj, so that Q m = ^/R s /Gj/2 = 10. In ac- 
cordance with above estimations, the widest synchroniza- 
tion domain spreads at the scale of ft itself. The widths of 
the domains decrease at much higher and much lower fre- 
quencies ojb, owing to decrease of K{ui). More details 
and finer steps can be seen in the right pane where the 
transresistance is plotted along the cut in Vf, — lb plane 
showing a typical devil's staircase curve. As a side note, 
the domains are not precisely single-connected, there is 
a fine structure of small "islands" of the same n,m near 
each domain. This structure is however too fine to be 
resolved at the scale of the plots. 

To address the quantum effects, we restrict ourselves 
to narrow synchronization domains where a new long 
time scale ~ (Aubj) -1 S> (u;^,/) -1 emerges. At this 
time scale, one can disregard the dispersion of quantum 




■Hi 



FIG. 3. Left: Synchronization domains (n/m) in the plane of 
normalized bias voltage and bias current. Right: Quantized 
plateaus of transresistance R = Vo/Io along the cut given 
by the line in the left figure. Dashed curve: continuous tran- 
sresistance as set by uncoupled Bloch and Josephson parts, 
R={-Kh/2e 2 )0jj{I b )/u B {V b ). 



noise and amplification coefficient and end up with a 
local-in-time action which is formally equivalent to that 
of a classical system subject to a white noise. Simi- 
lar approach has been applied to narrow Shapiro steps 
[19] , The slow variables in our case are the phases 
8(t),^f(t) of Bloch and Josephson oscillations, respec- 
tively. With those, the time-dependent current (volt- 
age) is represented as Io(t) = lo + Io(&Bt + 0(t)) 
(V (t) = V + I Q (ujt + *(*))), lo, V ■ We derive 
the effective action in the vicinity of the point in lb — 14 
plane where nuij — mujB = uj aiming to describing the 
(n, m) domain (In formulas for the action, h = 1 for com- 
pactness). 



S — Sb + Sj + S cp ; 

s B = r { dt (ee d - iT* B e' 2 d - (5uj B )e t 



(6) 
(7) 



Sj = g j dt[ ** rf - iT**l - (5uj)y d ) , (8) 

K\ 



S C p ~ u 2 



— J dt (- A B cos(m8 - + k) ■ 9 d (9) 
Ajcos(m9 - - k) ■ ^ d Y (10) 



Here, Sb,j describe Brownian motion of the phases 
in the absence of the coupling, g(r) » 1 being di- 
mensionless differential conductance (resistance), g = 
{h/Ae 2 )(dIb/dV ) (r = (e 2 /n 2 h)dV b /dI ), T} B - tiu 
being effective noise temperatures that depend on bias 
current and voltage. S sp gives energy (~ h | K |) gained by 
synchronization, n = axg(K). The coefficients Ab,j de- 
pend on lb, Vb as well as on n, m. We concentrate on the 
relevant variable 7 = m0 — n\I/ to reduce the action to 
the form 



S 



dt (a(77d - iT *l d - M - E c P sm77 d ) . (11) 
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FIG. 4. Left: Washboard potential for the phase difference 7. 
The hops over the barriers are synchronization errors. Right: 
The logarithm of the error rate across the synchronization 
domain. 



Here, the susceptibility a — gr/{gm 2 + rn 2 ), noise tem- 



perature T* = (T B m 2 g+T}n 2 r)/(gm 2 + rn 2 ), the energy 
barrier E cp — hw\ABnrK + AjmgK*\/(gm 2 +rn 2 ), and 
8u> = mSQB — nbujj. This action is formally equivalent 
to that of an overdamped particle moving in a trapping 
washboard potential V(pi) = — E cp COS7 — jhaSui (Fig. 4) 
and being subject to the thermal noise. If we neglect the 
noise, the motion obeys aj+dU( , Y)/d'y = 0. The station- 
ary solutions to this equation where 7 in trapped in one of 
the minima correspond to the synchronization of the os- 
cillations. They occur within a strip |&j|A<j = E cp /ha, 
in accordance with the estimations made. Beyond the 
strip, 7 increases with time corresponding to two unsyn- 
chronized frequencies. 

The synchronization errors are thermally-activated 
hops between the neighboring minima. Their rate de- 
termines the accuracy of the resistance quantization. It 
is clear that this rate is exponentially small, lnT = 
—E cp /T*, in the center of the synchronization domain, 
this guarantees the high quality of the resistance quanti- 
zation. The rate increases towards the strip edge owing 
to the lowering of the barrier in the washboard potential, 
lnT = -{E cp /T*)((l - y 2 ) + yarccos(y)), y = \6u\/Au. 
The coefficient E cp /T* ~ K depends on bias current and 
voltage as well as on n, m. We provide extensive illustra- 
tions of this dependence [21]. 

In FQHE, the excitation bear fractional charge/flux. 
The synchronization errors may also be considered as ex- 
citations at the background of a synchronization domain. 
One might conjecture that extra charge/flux induced by 
a hop over the barrier is fractional: this would be the case 
if the 2tt change in 7 is equally split between the phases 
6?, In fact, the situation is more complex since the hop 
takes a relatively long time ~ ah/E cp during which the 
charge and flux (related to the superconducting phase 
difference <p) may fluctuate. Owing to this, the average 
extra charge and flux transferred in course of a hop do 
not exhibit a strict quantization, 



However, in the limit g ^> r the extra charge approaches 
fractional value 2e/m, while the extra flux approaches 
— 1/n in the opposite limit. 

In conclusion, we have proven the feasibility of syn- 
chronization of Bloch and Josephson oscillations in super- 
conducting devices that is manifested as a (fractionally) 
quantized transrcsistance like in QH devices. High ampli- 
fication coefficient is required for the stable synchroniza- 
tion, this is readily achieved by using an LC-resonator 
with high quality factor Q. The minimum synchroniza- 
tion error rate is exponential in Q. 

We acknowledge fruitful discussions with J. E. Mooij, 
K. K. Likharev, and A. V. Zorin. Our research was sup- 
ported by the Dutch Science Foundation NWO/FOM. 
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Supplementary Material for 
"Quantum synchronization and transresistance quantization in superconducting 

devices" 

A. M. Hriscu, Yu. V. Nazarov 
Kavli Institute of NanoScience, Delft University of Technology, 2628 CJ Delft, The Netherlands 

Here we present supplementary material for our article " Quantum synchronization and transresistance quantization 
in superconducting devices". We give here details of the original Keldysh action that describes a general coupling 
circuit, specify classical saddle-point equations and details of their analysis. We present the derivation of simple 
action for narrow strips and complex formulas for the coefficients involved. We illustrate the dependences of the 
synchronization error rate for various synchronization domains on the parameters of the device proposed. 



DETAILS OF THE GENERAL ACTION 



The action for a general circuit where Bloch and Josephson junctions are present consists of a non-linear part 



S n i = I dt [ 2Ej sin (j) sin 



2hs sin q sin — 



the part describing the bias current and voltage 

^bias ~ 



- / dt 



eV b I b 

Id + 

7r 2e 



and the part describing the linear circuit. To get this part of the action , we use general expressions for Keldysh 
action of a linear system. Given a vector Xi of generalized coordinates, corresponding generalized forces ft, and the 
susceptibility matrix x{oo) that relates the coordinates and forces, 



the action can be written as (A = \ 1 ) : 

-/£(■ 



1 <? T a 



fXd 



The terms oc xxd describe the response as well as dynamics of x. The terms oc x 2 d describe the quantum and thermal 
noise. The quantity S is a correlator of random forces. By virtue of fluctuation-dissipation theorem it is expressed in 
terms of susceptibility as 



S u = l -(A- it) coth ( 



-) 

27V 



This matrix is positively defined (eventually, IraA is negatively defined at positive frequency) 

In our case, the generalized forces are the voltage source connected to the Bloch terminal and the current source 
connected to the Josephson terminal, such that the generalized coordinates are related to the charge and flux at the 
junctions. We use a common convention where positive currents are out of the terminals. We compute matrix U 
that relates sources and responses, those are current j in Bloch terminal and voltage v at Josephson terminal. The 
susceptibility matrix x — U /(—iiv). The explicit answer for U reads 



U 



V 

I 



u 



1 



R s Z 2 + (Rs + Z 2 )(Z 1 +Gj 1 ) 
Inverting the matrix, we obtain 

Z2 z± Z[ 



Z 2 + Zi+ Gj 1 
Z 2 Z\ 



-Z 2 Gj' 



-VjJ 



Rs 



Zl+Z2 Z2+Zl 
Zi t ^ _i 1 



Z 1 {R S Z 2 + Z 1 {R S + Z 2 )) 



R s + SR S K 
-K Gj + 6Gj 



(13) 
(14) 
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Here, like in the main text, we have introduced the amplification coefficient K = Z 2 /{Zi + Z 2 ) and coupling-induced 
modifications SR S = Z 2 Z 1 /(Z 2 + Z x ), 5Gj = l/(Z 2 + Z{). 

The variables <p,q differ by factors from the generalized displacements $ and Q, <f> = 2e$, Q — (e/w)q. With this, 
the part of the action representing the linear circuit reads: 

S = 5i + S 2 + S 3 + S 4 (15) 

(«*-«.-*"«>•) (») 

f duj cu , / ui \ Imif(w) / j , j , N 
54 = 7 2^ 2 C ° th (2?) * 2tt ^""^ ~ (19) 

with G = Gj + 5Gj, R = R s + 8R S . 

CLASSICAL EQUATIONS 

The classical equations for the circuit can be obtained either by varying the action with respect to (f>d,qd and 
setting (j)d 7 qd = or by applying Kirchoff equations to the circuit. The latter will allow us to keep the equations in 
the differential form 

The classical equations for q and <p ar e separate in the absence of the coupling, 

/ + / c sin^ +2 A-f .0 (20) 

V + V c smq+-R s ^ =0 (21) 
7r at 

and are obtained by applying current conservation in Joscphson branch and summing up the voltage drops along the 
Bloch branch. With the coupling elements, these equations include the current in the capacitor I cap and voltage drop 
Vl over the inductor, 



2eRj dt 



/ + 7 c sin</»+^^^-/ cap = (22) 



V + V c smq+ e R s ^+V L =0 (23) 
tt at 

that are related to the corresponding current and voltage by I cap — CV, Vl = LIl- Two extra equations are obtained 
by applying current conservation in the node connected to the capacitor and the inductance, and summing up the 
voltage drops across the Josephson junction, capacitor, and inductance, 

^4 - leap +I L \ ~ Vca P + Vl - ( 24 ) 

We reduce these 6 equations to 4 equations that express time derivatives of q, 0, Vl, II in terms of these 4 variables. 
We solve these evolution equations numerically by a rk4 solver. A typical computer run is as follows. We fix the 
device parameters Ej, Eg, Rs,Gj, L,C for the whole run. We vary changeable parameters h,Vb with small steps, 
cither along a line in the I), — Vt, plane or scanning a square in the plane. For each point in the It, — Vj, space we start 
with an arbitrary initial condition and "wait" (typically, 30 — 40 Q m /fi) for equilibration of the resulting oscillations. 
Then we check whether the resulting orbit is a periodic one in the space of q, <fi, Vl, II that is, if initially at the point 
{q, <f>, Vl, II), it arrives to the point (q + 2irn, <fi + 2i\m, Vl, II) with some integer n, m after a time interval. If the 
periodicity is found, the program outputs the numbers n and m and the point in the lb — V, plane is marked as to 
belong to the synchronization domain n, m. Otherwise, the point is regarded as to belong to the chaotic domain. In 
a faster version of the simulation, Ib,V> are continuously and very slowly updated with time, and the periodicity of 
the orbits is monitored constantly. 
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In this way, we have obtained the plots given in the main text. Those are made with relative resolution 10~ 4 in 
lb, Vb and show single-connected synchronization domains at the chaotic background with widths quickly decreasing 
with increasing max(n, m). 

If we concentrate at a domain boundary, say, of the (1,1) domain and increase resolution, we find more structure. We 
are able to see the "paddles" of the chaotic domain within (1, 1) domain, the small "islands" of (1, 1) synchronization 
in the chaotic region and even smaller synchronization domains with large n, m. The structure exhibit fractal self- 
similarity upon scaling the resolution and the size of the region scanned. This is what is generally expected from the 
transition between commcnsurability and chaos and should not surprise. Since the structure is seen at high relative 
resolutions only, we expect it , in distinction from wider domains, to vanish at any realistic noise level and therefore 
did not investigate it in detail. 



DERIVATION OF THE SIMPLIFIED ACTION 



Our goal is to derive a simplified action for slow variables that are the phases <!/, 9 of the Josephson and Bloch 
oscillations, correspondingly. One can draw a similarity with a well-known phenomenon of Shapiro steps in Josephson 
junction whereby Josephson oscillations are synchronized with an external a.c. current signal. Our case can be 
regarded as a sort of spontaneous emergence of such signal whereby the synchronizing signal for Josephson part is 
produced by the Bloch part, and vice versa. To this end, we start with deriving the action for Shapiro steps in the 
resistively shunted Josephson junction. 



Simplified action for Shapiro steps 

Our starting point is the full action for the resistively shunted Josephson junction, 



= J dtl 2Ej sin sin ^ - 1 + ^ 4>d + M>d^\ + S n (25) 



where the noise term S n we write at the moment as 

s n = -i\ J dtdt'MMt - t')Mt') (26) 

without specifying the kernel s. 

We assume G ^> e 2 /fi, (effective) noise temperature <C Ej, and small a.c. amplitude I <C Ic- Under these 
conditions, the realizations of <f>, 4>d are close to solutions of the stationary saddle-point equation, that is, one without 
noise, and cf) d is typically small. The specifics of the situation that these solutions are time-dependent periodic 
oscillations and therefore are not unique: For each solution <j>(t) there is a phase-shifted solution <j)(t + ^f/uj), oj being 
the oscillation frequency. The field '3/ can therefore be regarded as a Goldstone mode. To come to the effective action 
for one allows for slow time-dependence of this field, expresses <j>(t) , <j>d(t) in terms of ^(t),^d(t), and substitutes 
these expressions into the action. This rather straightforward program has some less trivial implementation details 
outlined below. 

The classical equation without noise is obtained by varying the general action with respect to (f>d- It reads 

£jSin</> -^ + ^ = (27) 

To solve it, we first make it dimensionless by introducing Ic = 2eEj, dimcnsionless current j = I/Ic, and dimcn- 
sionless time r such that 

dr Ae 2 E.,dt V ; 

With this, the equation reads 



4> + sin 4> — j = 



(29) 
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To solve it, let us change the variable to Z = e 1 ^, Z — iiftZ, so the equation becomes 

Z+ \ (* - 1 - 2ijZ) = - z2 _f_ 2ijZ = -dr/2 (30) 

Let us substitute j = \{y + y~ r ), thus defining a convenient parameter y > 1. With this, we can factorize the 
denominator, 

Z 2 - 1 - 2ijZ = (Z- iy) (Z - ly- 1 ) (31) 

Integrating both parts we get 

In ( Z ~ iy \ - 1 - = -r/2 + C (32) 

\Z-iy~ 1 ) i{y-y- 1 ) ' v ; 

We choose the constant C in such a way that Z = 1 at r = and exponentiate to obtain 

Z — iy 1 — iy -1 



Z — 1 1 — iy 



exp (— irujj) (33) 



with the dimensionless Josephson frequency vj — \{y — y 1 ) = \J j 2 — 1. The frequency of correct dimension is then 

We can now rewrite the equation for Z as 

Z — iy t 

- ^— = y exp(-irwj - «x) = V x (34) 

Z-iy 1 



introducing a phase factor 



to obtain 



s + 

cxp(ix) = -i- - 

1 - zy 



X - y . 1 - y 2 i . 2oj. 



= * — — "T = W + i ', = i 



yX-l y y-l/X y yX - 1 
That we can expand in 1/X, that is, in harmonics: 



i \ - 1 

Z = - - i2uj > — 

" n=l y 



There are only harmonics with negative and zero frequency. 
Let us compute the voltage dcfi/dr. It is given by 

Z ■ dZ _! <9Z _! 2 yX 2 y 

■> (x y)( y x 2vj v 2 + 1 - yjx + i/x) 

It is a less expected property of this equation that the inverse voltage has only zero and first harmonics not involving 
any higher Fourier components. 

(d^/dr)- 1 = cuf (j - cos( Wj t + X )) • (35) 

Derivation of slow-variable action 

Let us use the solutions found. We can express the classical field </>(<) as 

0(t) = F(t + *(t)/wj). (36) 

F(t) being the solution found above. It is important to notice that a change of the field ^ amounts to a slow part of 
the change of the superconducting phase <f>, since the change of <j) per period of Josephson oscillations is 2n. 
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To write a similar expression for the quantum field 4>d is less trivial. One might conjecture that Eq. 36 is valid for 
both parts of the Keldysh contour, that is, for <j> (t), and get the expression for <fid in this way. This conjecture is 
however wrong not satisfying the saddle-point equations. To find the true dependence of 4>d at the time scale of the 
oscillation period, we need to inspect the "quantum" saddle-point equation obtained by varying the full action with 
respect to <f>. Assuming small 4>d, the equation in terms of dimensionless time r reads 



To solve this equation, we notice that 



From the other hand. 



i cos 6 — 6d = 



^- \n(<j) d ) = cob(0). 
dr 



(j) + sin <f> — j = — » <f> + <j> cos (j> — > cos (j) = — — In (j)^J (37) 

We conclude that 4>d = C '(dip / 'dr)^ 1 . We fix the normalization constant C from the condition that the slow variations 
of 4>d correspond to variations of ^d, so that the average of <fid over the oscillation period cf>d = ^d- Using the explicit 
formula for difi/dr, we obtain 

Mt) = *d(l - d~ 1 ™$(u>jt + X + f)) = y d (d^/dry 1 -4 (38) 

J 

Now it is time to substitute (j>(t), 4>d to the action and average it over the period. In zeroth order, this cancels the 
action: indeed, this is generally expected from Goldstone fields. We deal with the residual terms one-by-one. The 
conductance term where we differentiate VP when taking the derivative with respect to time, reads as follows 

where G* — Gvj/j is, logically enough, is the differential conductance (dV/dlb)^ 1 . Next we look at the terms that 
come with the external current. We expand the current in harmonics of frequency w k wj, 



/(*)=/+ 5^ Re (J n e inut ) 



With this, the action should become 



(A2) + ±£ Re (/„(&);)) 

where AI = I — I(u>) is the external current minus the current corresponding to the frequency w and (4>d)n are 
the normalized harmonics of <fid- It is a rather peculiar property of the model in use that (4>d)n — except n = 1, 
= ~j exp(i(^ + x)). Owing to this, the model predicts Shapiro steps only if the external frequency equals to 
Josephson frequency. 

By no means this is a general situation. Generally, one expects the Shapiro steps at all multiples of Josephson 
frequency. In fact, this drawback of purely resistive shunting is well-known in the field of Josephson dynamics and 
arises from the absence of frequency dependence in Gj. If a shunting capacitor is added to the model to provide the 
frequency dependence, such Shapiro steps readily emerge. We see that also in our simulations, where the frequency 
dependence is provided by the coupling LC circuit. To simulate the generic situation without a complication of the 
model in use, we will force the presence of actual harmonics in <pd at all frequencies 

(4>d) n = -a n r n exp(m(* + x))- (40) 

The values of a n ~ 1 at n ^ lare subject of choice and can be regarded as extra model parameters. 
The noise part of the action after the substitution takes the form 




lY, S (™MMn\ 2 )- (41) 
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So we need the noise kernel at multiples of the Josephson frequency. In accordance with the fluctuation-dissipation 
theorem, s(uj) — ^wcoth(w/2T). We assume quantum limit wj » T, so that S(uj,j) = j^\ujj\. 

We can rewrite this contribution to the action as ^f^2T*, introducing the effective temperature of the Josephson 
oscillations, 

T * _ U J 
J ~ Wf=T' 

Collecting all terms, we obtain the effective action in the following form 



f , / 9U T • G* nrivt G* * 2 

J '"{ ' ,ln|; "'- '' 2r < - ) (42) 



l 4e 2 J 4e 2 2 

where the effective potential U(^) gives the interaction with the a.c. current, 



dU 

9* 2e 



1 ^(A7) + l - J2 ^ (/^nj-" exp(-in(* + X )))^ . 



In distinction from the original action, this action is local in time. It can be solved by integration over subsequent 
time slices. The resulting evolution equation is a Fokker-Planck equation for distribution of P(4 r ). For us, it is 
sufficient to note that the action is equivalent to that of an overdamped classical particle subject to white noise with 
effective temperature T* . The synchronization errors correspond to thermally activated hops over the barriers in the 
potential {/($) and the estimation of their rate is given by Boltzmann factor exp(— AU/T*). 

Action for the Bloch part 

We proceed in a similar way to obtain the action for the slow variables 9, 9 d of the Bloch oscillations. 
The starting point is the action in terms of Q,Q d 



2E S sin Q sin 9±-±V(t)Q d + QQ^-^- ) + S n (43) 

A IT : ~ 



with the noise term that we write at the moment as 

S n = -i^J dtdt'Q d (t)s(t - t')Q d {t') (44) 
The classical part of Q is expressed as 

Q(t) = F(t + 0(t)/u B ), 

where ujb is now the Bloch frequency. The slow field 9 that adds to the phase of the oscillations. The quantum part 
is obtained as 

Qd(r) = 9 d (l - 1 cosset + X + 9 d )) = 9 d (dQ/dr)- 1 ^ 

where we introduce dimensionless bias voltage v = Vb/Vc and dimensionlcss Bloch frequency vb = V« 2 — 1- 

Again we need to substitute the above expressions Q(t),Q d to the action and average it over the period. We 
proceed as in the case of Josephson junction. Let us explicitly concentrate on the dissipative part. According to 
fluctuation-dissipation theorem, S(uj) = ^^ L ujcoth(uj/2T). Assuming ujb ^ T, we can rewrite the action as 

(ft eft eft c>* 

Vd S * = v_d 2T *^s (45) 

2 2 7T 2 K J 

introducing effective temperature 

UJ B 



T* B = 



2 
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With this, the effective action reads 



s= l<i'l-§o d + ee^-i2T* e ^kk) 



where the effective potential U (9) describes the synchronization with the a.c. voltage signal, 

\ (^AV) + ^J2 Re ( V ^- n exp(-in(9 + x)))^ ■ 



dU 
~d9 



Here, we add /3 n that are the harmonics of Qd at n ^ 1 describe the synchronization at multiples of Bloch frequency. 
The differential conductance and effective temperature read 

H s = K s — ) 1 =- — / „ ■ 
v 4vVv 2 - 1 



Coupling the parts 



Now we can combine the previously derived slow- variable actions for Bloch and Josephson p arts . The only remaining 
work to substitute the harmonics of Qd, <f>d, Q, into the coupling part of the full action (Eq. 17). We concentrate on 
the vicinity of a curve in the lb — Vb plane where nujj(Ib) = muiB{Vb) = w. In the vicinity, we can disregard the details 
of frequency dependence of K(uj) replacing it with a (big) complex number K. In principle, the synchronization can 
be achieved at all multiples of the frequency uj. This would however involve K{uj) at multiples of a;. Since we assume 
that K is big only in the vicinity of the resonant frequency, we can safely disregard these terms. 

The harmonics we need are 



(f) d {t) = - — ^ d gi«("jt+*+X./) _|_ p -in{uJ.jt+^+x.j) 



Qd(t) 



' 2v r 



e im(ui B t+8+XB) _j_ e -im{u) B t+8+XB) 



(1) _ U J ( e in(u,jt+<Sr+ X j) i e -in(wjt+V+Xj)'\ 
{] + VJY \ ' > ' 

(Q\ = ( e im(u,Bt+$+XB) + e -irn{u B t+e+ X B) 

(v + V B ) n \ 



(47) 

(48) 
(49) 

(50) 



With this, we arrive at 



S cp — / dt(-A B cos(m9 - + n)9 d + Aj cos(m9 - - n)^ d ) (51) 

2ll J 

with 

J mj n (v + Vb)™' B nv m (j + vj) n ^ 

The phases xj,b are irrelevant and can be canceled by corresponding shifts of , 5, 9. This does not apply to the phase 
k of the amplification coefficient. 

When deriving the action, we have disregarded £4 that describes the correlation of the noises induced by coupling. 
This term is small in comparison with the noise produced by Rs, Gj provided SR$, SGj <C Rs, Gj: this is what we 
assume. 
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Reduction to the single variable 

The resulting action 

S = S B +Sj+S cp ; (53) 

s B = rjdt (ee d - iT* B e 2 d - (5u B )e d ) , (54) 



S B=gJ dt - iT*^> 2 d - {5ujj)^ d ) , 



(55) 



s -J Kl 



J dt(-A B cos(m6 - + n)6 d + Aj cos(m8 - - k)^). (56) 



depends on two pairs of variables while the only variable that enters the action in a non- linear fashion is 7 = mO — n^. 
To reduce the action to that of this relevant variable, we substitute 

(57) 

* d = -^--; e=f---; (58) 

2n n 2m m 

introducing the fields 7,7,2, &,&d- The resulting action is quadratic in <r, cr<2 so that this field can be integrated out. 
Since the integral is Gaussian, the integration is equivalent to finding an optimum in a, a d and substituting the optimal 
values a, a d found back into the action. 
In this way, we arrive at 



7 


a 

--; < 


? = +^- 


a 


2n 


n 


2m 


m 


Id 


a <i. 


e = ^L- 


<?d 


2n 


n 


2m 


m 



S = j dt (a(77d - «T*7^ - 5w) - E cp smj-j d ) . 



(59) 



Here, the susceptibility a = gr/(gm 2 + rn 2 ), noise temperature T* = (T B m 2 g + T}n 2 r)/(gm 2 + rn 2 ), the energy 
barrier E cp = hw\A B nrK + AjmgK*\/ (gm 2 + rn 2 ), and Sco — mSui B — nSuij. Expressing changes of the fields 9 in 
terms of 7 we arrive at the relations 

Sq _ mg <^> _ -nr 

2e gm 2 + rn 2 ' 2ir gm 2 + rn 2 

that give the change of the charge and flux passed in the output circuits upon hopping over the potential barrier. 



NUMERICAL ILLUSTRATIONS 

In the main text, we have argued that high quality factors Q, or, equivalently, large amplification coefficients K are 
required for exponentially good synchronization, lnT oc \K\ < Q. The effective quality factor is limited by y/RsGj/2 
and can therefore be sufficiently large under assumption of well-defined oscillations Rs 3> fi/e 2 , Gj » e 2 /h. It is 
important to find numerical estimates on proportionality coefficient in the above relation. Apparently, the coefficient 
depends very much on the parameters of the model and tend to be very small for exotic synchronization plateaus 
with high n, m. Here, we present these numerical estimates. 

We concentrate on evaluating a dimensionless coefficient A defined by 

r.exp(-^ 

T estimating the synchronization error rate in the middle of a synchronization domain. Thus defined A does not 
depend on the characteristics of LC resonator. It follows from the formulas given that 

^ _ E cp 2ir _ yj ' A 2 B n 2 r 2 + A 2 B n 2 r 2 + 2AjA B nmgr cos(2k) 
T* \K\ T_B_ m 2g _|_ ^n 2 r 
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X being the phase of the amplification coefficient. All coefficients Ab,j,T b j,r,g depend on n,m and dimensionless 
voltage and current biases j = h/I c -> v = Vb/Vc, 

r; ^ = —— -%==—-; (62) 



mj n (v + Vv 2 - l) m ' nv rn (j + yjj 2 - l) r 

r = f Vv2 ~ 1 ; g = g ^fE± - whevef=(e 2 /n 2 h)R s ;9 = (h/Ae 2 )Gj (63) 
v j 

T A = 1 . Tl = 1 . (64 ) 

huj Amv^Jv 2 - 1 ' hw Anj^Jj 2 - 1 

We introduce two dimensionless parameters that characterize asymmetry of the Bloch and Josephson parts, 

X = f/g,Y = Ej/E s . 

The symmetry is achieved at X = Y = 1. 

The biases are related by synchronization condition ntuj = muiB that in dimensionless notations reads 

y/v*-l = XY-y/p^l. 

m ' 



Main synchronization domain n = m = 1 



Let us concentrate first on the main synchronization domain n = m = 1. We plot in Fig. 1 the coefficient A versus 
the dimensionless bias current j. 

Upper left plot shows the coefficient for different cos(2k) = 1,0.5,0,-0.5,-1.0 under symmetry conditions X = 
Y = 1. At any j, the coefficient reaches the maximum for purely real K(cos(2n) = 1) and vanishes for purely 
imaginary K, this being an artifact of the symmetry of the setup. For all cos(2k) A is zero at the threshold j = 1 and 
monotonically increases saturating at big j. This is explained by a rather peculiar dependence of T* /ftw on the bias. 
This quantity diverges near the threshold and quickly drops with increasing the bias. This drop compensates decrease 
of E cp at big biases. We note that the strip width that does not depend on T* becomes smaller upon increasing j. 

Upper left and lower right plots show the influence of asymmetry on the coefficient, for different Y and X, respec- 
tively. In all cases, asymmetry decreases the coefficient, and A vanishes in the limit of either big or small X as well 
as in the limit of either small or big Y. The reason for this is that at big asymmetries the synchronization condition 
implies that one of the parts is biased close to the threshold where the effective temperature is high. 

Asymmetry improves A in the case of purely imaginary impedance when it vanishes in the symmetric case. This 
is illustrated in the lower right plot where it is shown that A attains a relatively big value ~ 0.5 for not so big 
asymmetries. 

The maximum A = 2 is achieved in the case of real impedance, symmetric conditions and large bias. 



Domains n = l, m/1 or m = l, ra/1 



Let us consider synchronization domains with higher n, m. Among those, the case with either n = 1 or m = 1 is 
rather special. In this case, we also expect steps in our high-damping limit model without adjusting a, ft, so we set 
a p = (3 p — except p = 1. In this case, either Aj or As equals zero, and the coefficient A does not depend on the 
phase of the amplification coefficient. We plot in Fig. 2 A versus j for a variety of n, m (note log scale of the vertical 
axis). We consider only symmetric setup X = Y = 1. The behavior of the coefficient in n = l,m and m = l,n 
domains is rather similar so for concreteness we concentrate on m = l,n. As in the previous case, the coefficient 
vanishes at the threshold. Unlike the previous case, it drops upon increasing the bias current. The reason for this 
is that E cp is proportional to the n-th harmonics of the oscillation that scales as j n at large bias. The coefficient A 
therefore reaches the maximum at some j n . With increasing n j n approaches the threshold, while the value of the 
maximum quickly drops. Roughly, it is decreased by a factor of 3 if n is increased by 1. 

The maximum A found is w 0.1 for n = 1, m = 2 and m = 1, n = 2 domains. 



14 




1 1.5 2 2.5« 3 3.5 4 1 1.5 2 2.5 3 3.5 4 



3 j 

FIG. 5. The coefficient A for main synchronization domain n = m = 1 versus the dimensionless current bias j. Upper left: 
symmetric setup X = Y = 1. The phase of the amplification coefficient takes values cos(2/-c) = 1,0.5,0,-0.5,-1.0 from the 
upper to the lower curve. Upper right: real K, X — 1 and asymmetry parameter Y takes values 1, 2, 0.5, 3, 0.33 from upper to 
lower curve at large j. Lower left: real K, Y = 1 and asymmetry parameter X takes values 1, 2, 0.5, 3, 0.33 from upper to lower 
curve at large j. Lower right: Purely imaginary K , Y = 1, X takes values 0.5, 2, 0.33, 3 from upper to lower curve at large bias. 

Domains n / l,m ^ 1 

For completeness, we illustrate A for two general domains: n/m = 2/5 and n/m = 3/4. We set a = ft = 1 for all 
harmonics. 

The coefficient A approaches zero at the threshold and peaks rather close to it dropping down at larger bias currents. 
The dependence on the phase of the amplification coefficient is bigger for n/m — 3/4 step since this ratio is close to 
1 where in symmetric case the cancellation of A takes place at purely imaginary amplification coefficient. In both 
domains, maximum value of A ~ 0.15. 
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FIG. 6. The coefficient A for synchronization domains n = l,m/ 1 and m = l,m/l versus the dimensionless current bias j. 
The coefficient does not depend on the phase of the amplification coefficient, and we assume symmetric case X = Y = 1. Left: 
m = 1, n changes from 2 to 6 from upper to lower curve. Right: n — 1, m changes from 2 to 6 from upper to lower curve. 




FIG. 7. The coefficient A for synchronization domains n/m = 2/5 (left) and n/m — 3/4 (right) versus the dimensionless 
current bias j. We assume symmetric case X = Y = 1. cos(2«:) takes values 1, 0.5, 0, —0.5, —1 from upper to lower curve. 



